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I It has been observed that the statistical distribution of the eigenvalues of ran- 

dom matrices possesses universal properties, independent of the probability law of the 
^ • stochastic matrix. In this article we find the correlation functions of this distribution 

in two classes of random hermitian matrix models: the one-matrix model, and the 
. two-matrix model, although it seems that the methods and conclusions presented here 

! will allow generalization to other multi-matrix models such as the chain of matrices, or 

^ I the 0{n) model. We recover the universality of the two point function in two regimes: 

short distance regime when the two eigenvalues are separated by a small number of 
other eigenvalues, and on the other hand the long range regime, when the two eigen- 
values are far away in the spectrum, in this regime we have to smooth the short scale 
5^ I oscillations. We also discuss the universality properties of more than two eigenvalues 

correlation functions. 
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It has been observed experimentally [|I], for disordered systems that the distri- 
bution of the energy levels of such systems is universal in some regime. For instance 
the connected correlation function between two levels separated by a small number 
of other levels does not depend on the system, while the density of levels is very de- 
pendent of the specific details of the system. Such a phenomenon is also observed in 
numerical simulations of large random matrices |^ . The consequence of this is that the 
correlation functions can be obtained from a gaussian model. This conjecture has been 
proved in the special case of a one hermitian matrix model 0, |§] . It has also been noted 
that the connected correlation functions of more than two eigenvalues should present 
also a stronger universality than the density itself. Brezin and Zee Q have calculated 
explicitly correlation functions of eigenvalues of a class of stochastic hermitian matrices 
of large size N. They have found that some statistical properties of the eigenvalues are 
universal in the large limit, and can thus be obtained from the correlation functions 
of the gaussian model. In addition, they have discovered that the two-point correlation 
function, after smoothing of the short scale oscillations, is universal while all other 
smoothed correlations vanish at the same order. This property could be understood 
by a renormalization group analysis |^ which indicates that the gaussian model is a 
stable fixed point in the large N limit. 

Random matrix models appear also in another field of physics: they have been 
discovered to describe a regularized theory of random surfaces and string theory [|, ^ . 
The string theoricians have thus developed very powerful tools for the study of the 
critical properties of these models, as orthogonal polynomial methods, saddle point 
approximations or loop equations. 

The analysis of |Q is based on the, by now standard method of orthogonal polyno- 
mials. An essential ingredient in the final answer is a proposed ansatz for an asymptotic 
form of the orthogonal polynomials P„ in the limit N oo and N — n finite. In the 
ansatz is verified in the case of even integrands, and only up to an unknown function. 
They also examine the smoothed correlation functions in the large scale regime, using 
a diagrammatic method. 

Here, we propose a direct proof of the ansatz for the orthogonal polynomial, using 
a saddle point method, which does not depend on the parity of the integrand, and 
which allows to determine the previously unknown function. Moreover, we extend 
this method to another class of matrix-models known in 2D gravity as the two-matrix 
model. This model presents a larger class of critical points than the one-matrix-model, 
it is therefore interresting to find if it presents the same universality. 

In addition, using the loop equations approach [1^, |ll|, we can generalize the di- 
agrammatic method of |^ to find all the connected correlation functions of the one- 
matrix-model, and discuss their universality properties. This recursive method allows 
explicit calculations, and we write the three and four point function. This method 
could also be generalized to another class of matrix models: the 0{n) model, which 
has the same kind of loop equations and topologial expansion procedure ||T6[| . 
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1 The 1-Matrix-model 



Let us first explain the problem and recall the method used in ref. to explicitly 
evaluate the eigenvalue correlation functions. 

We consider a N x N stochastic hermitian matrix M with a probability law of the 
form: 

, , 1 -A^trV^(M) , ^ 

V{M) = -e ^ ^ (1.1) 

where V{M) is a polynomial, and Z the normalization (i.e. the partition function). 
We want to derive the asymptotic form for large of various eigenvalue correlation 
functions. All of them can be obtained from the correlation functions of the operator 
(9(A): 

1 1 ^ 

0(A) = _ tr 5(A - M) = - ^ 5 (A - ^.) ( 1.2) 

(the fii being the eigenvalues of M). Indeed, for any set of functions (/i, ...,fk), one 
has: 

( tr /i(M) . . . tr MM)) = J dXJ,{X,) . . . dX,/fe(A,) (0(Ai) . . . O(Afc)) . 

The correlation functions of the operator 0{X) can in turn be expressed in terms of 
the partially integrated eigenvalue distributions like p(A) the density of eigenvalues, 
which is the probability that A belongs to the spectrum of M, p2{X, fi) the probability 
that A and fi are simultaneously eigenvalues of M and more generally p„(Ai, . . . , A„) 
the probability that Ai ...A^ are simultaneously eigenvalues of M. For Ai 7^ A2 . . . 7^ A„ 
(else some additional contact terms have to be added) we find 

{0{X,)0{X,) . . . 0(A„)) = -^j^;f^Pn{Xu • • • , A„). 

Actually the interesting functions are not directly the 0{X) correlation functions, 
but their connected parts. Indeed, at leading order, when A^ ^ 00, we have the 
factorization property 

(0(Ai)0(A2)...0(A.))~ n (^(^0). 

i=l,...,n 

and thus no new information can be obtained from the complete n-point function. The 
connected function which will be denoted 

(0(Ai)...0(A„)U, = a(Ai,...,A„) 

The method of orthogonal polynomials allows to determine directly all these connected 
functions from only one auxiliary kernel k{X, fi) = K,{p, X) (see refs. |]T^ ^ [T^ or 
appendix 2 for details). For instance we have: 

Oi(A) =p{X) = KiX,X), 

02iX,p) =-K\X,fi), 

OsiXu A2, A3) = 2/t(Ai, X2)k{X2, Xs)k{Xs, Ai), 
Oi{Xi, A2, A3, A4) = -2k(Ai, A2)k(A2, A3)k(A3, A4)k(A4, Ai) 
-2k(Ai, A3)k(A3, A2)k(A2, A4)k(A4, Ai) 
-2k(Ai, A2)k(A2, A4)k(A4, A3)k(A3, Ai), 
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and analogous expressions for larger values of n: 



a = (-ir^ E 11'^ (^-^^J (1-3) 

cyclic permutations cr i=l 

The kernel k is simply related to the polynomials Pn orthogonal with respect to the 

-NV{X) 
measure dA e : 

1 N-l 

'^(^'/^) = 77 E Pn(A)P.(/x)exp[-(iV/2)(\/(A) + Vifi))] ( 1.4) 

n=0 

With the help of the Darboux-Christoffel theorem (Appendix 3) we can rewrite it: 

.(A,/.) oc - PN-.WPNif.) ,,p[_(^/2)(V(A) + Vim- ( 1.5) 

So, the asymptotic evaluation of correlation functions is reduced to an evaluation of 
the kernel k(A,/x) and thus of the orthogonal polynomials Pn(A) in the peculiar limit: 
large, N — n finite, and A G [a, b] ( [a, b] being the support of p(A) ). 
The ansatz proposed in @| in the case of even potentials V{M) (for which b = —a) 
was: 

P„(A) (X e^^^^^/^^^cos(ArC(A) + (iV-n)¥.(A) + x(A)), 



where 

A = a cos (p , 
/(A) = asiny9 , 
^C(A) =-vrp(A), 

x(A) remaining undetermined, except in the gaussian and quartic cases, for which: 
X(A) = - 7r/4. 

We will prove below that this ansatz is still true for any V, and that x is always 
of the form x = y)/2 + const. The method which leads to the proof is also interesting 
in itself because it uses some general tools of the saddle point calculations of the one- 
matrix model 0. This method will be extended to the multi-matrix model in the 
next section. 

1.1 Orthogonal Polynomials 

We thus have to consider the set of polynomials P„ (n is the degree), orthogonal with 
respect to the following scalar product: 

r -NV(X) 

(P„ .P^) = 5nm = j P„(A)P™(A). ( 1.6) 

Remarkably enough an explicit expression of these orthogonal polynomials in terms of 
a hermitian matrix integral can be derived (see appendix 1, or p5|): 



P„(A) . /d-Me-"^" det(A - M), ( 1.7) 
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where M is here a n x n hermitian matrix, and Z„ the partition function ( | corre- 
sponding to n X n matrices. We will use this expression to evaluate P„ in the relevant 
limit, i.e ^ 1 and N — n = 0(1) by the steepest descent method. Let us first 
consider the matrix integral: 



Z{g,h,X) = J (T'^'^M e 



[n/g) tr V(M) 



V{z) =V{z) -h\Ti{\- z). 

With these definitions Pn(A) is proportional to Z{g = n/N, h = 1/N, A), and thus we 
need Z or equivalently the free energy Fn = —n^"^ In Z for h and g — 1 small. 

As we are interested only in the A dependance of F, let us differentiate F with 
respect to A: 

— uj[z = X,g,h,X) 



d\ g 



where uj(z) is the resolvent: 



u;{z) = -(tT^—). (1.8) 
^ ^ n\ z-Ml ^ ' 

Since we want the asymptotic expression of uj{X)n'^h/ g, we need uj up to the order 
1/n. It is known from the random matrix theory , that F„, and also uj have an 
expansion in powers of 1/n^ which is the topological expansion. At order 1/n, only the 
contribution of the sphere is required. We thus replace uj by its dominant contribution 
obtained by the saddle point method. With this approximation, uj may be written: 



n . 1 z — Xi 

j=l,...,n ' 



where the Aj are the eigenvalues verifying the saddle point equation: 

V(A.)4i: ' 



g rij^.\i- Xj 

i.e. they extremize the integrand of 

NViXi 



Z{g,h,X)^ J fie ^ ''dXi Ui^^-X 



\2 

i=l i<j 



In the large n limit the Aj are distributed along an interval [a, b] with a continuous 
density p(A) [|, g. Then: 

uj{z)= [\p{^)^ (1.9) 

Ja Z — /i 

and the saddle point equation becomes: 

^V'(/i) = w(/i + zO) +cj(/i-iO) for /iG[a,6] (1.10) 
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Note an important property of this equation: At a and h fixed it is linear and tlierefore 
the derivatives of uj with respect to g or h will also satisfy a linear equation. 

At leading order we introduce the resolvent uJo{z) = uj{z, g = 1, h = 0, X), and write: 

gu{z, g, h, A) = uJo{z) + {g - l)^g{z, A) + hQh{z, A) + 0{l/n^), 

where we have defined the two functions: 



dguj 



dg 



dguj 



g = l,h=0 



dh 



3=1, h=0 



The function ^0(2;) is the resolvent of the usual one-matrix model, and from eqs. 
1.9|) and ( | 1.10|) we obtain: 



ujq{z ± iO) = ^V'{z) T iT^p{.z) 



for z G [a, h]. 



As we noted above, the two functions and Vt^ obey linear equations, obtained by 
differentiation of the equation ( | 1.10|) satisfied by uj. Following a method introduced in 
ref. |jl5| we can then easily determine them from their analyticity properties, and the 
boundary conditions. VLg[z) verifies the linear equation: 

VLy{z + + ng{z - iQ) =0 



and behaves as I/2; when z ^ 00, and as 1/ \/z — a\J z — b near the cut end-points a 
and b, because u behaves as z — aV z — b. These conditions determine ilg uniquely: 

ng{z, A) = . ^ ■ 

^{z-a){z-b) 
The same method applies to VLh which satisfies 

Qh{z + iO) + nh{z - zO) = 1/(A - z), 



and behaves like Vth ~ Oil/z^) for z large, VLh ~ — a^/ z — b near a, 6, and is 

regular near z = \. It follows: 



Vth{z,\) 



1 



1 



V(^-a)(z-6) 
Then if we set 2; = A: 

nh{z = A) 



^{z-a){z-b)-^{\-a){\-b) \ 



z-X 



^ ln^(A-a)(6- A). 



2^(A-a)(A-6) 2 dX 
We now have the necessary ingredients to determine dF/dX: 
BF 

-n"— = Ntuoiz = A) + (n - N)ng{z = A) + n^iz = A) + Oil/N). 
oX 



We still have to integrate all these terms with respect to A. In order to integrate 
^0 = y 1'^ — iT^p, we introduce a primitive of p(A): 



C(A) = -7r /'dX'p(A'). 

J a 



We also need a primitive of VLg. For this purpose, we parametrize \ = \{a + h) — \{h — a) cos(^, 
so that VLg (xl/ sin ip and its primitive is simply 

11 

/ f2g(A')d\ , = V? where A = -(a + 6) — -(6 — a) cos(/9 

J a Z Z 

Finally, the result takes the form: 

-n^F{\±iQ) = I W(A ) ± iNCjX) t i {N - n)ip ± 
— \ In ^JJ}^^-a){h^^X) + const . 

Since P„ is a polynomial, we have P„(A) = | [Pri(A + iO) + P„(A — iO)] (this comes 
from the fact that there are actually two saddle points giving contributions of the same 
importance), and therefore: 



^ NViX)/2 1 

"n(A) = \ -e I cos 

./7(A) 



N({\) - {N - n)ip{X) + i(y9(A) + const 



1.11) 



where /(A) = y^(A — a)(6 — A). The constant factor \j2/n is fixed by the condition 
that k(A, A) = p(A). In the case of even potentials V, parity considerations imply that 
the arbitrary constant phase is const = —n/A. For general potentials the constant 
phase remains undetermined at this order, but we note that the general form of P„ 
does not depend on the parity properties of V. 
We have thus proven the ansatz needed by p| . 



1.2 Connected Correlation Functions 

From this asymptotic expansion of P„ one can now derive the kernel ( | 1.5| ) k{X, /i), and 
then the connected correlation functions, in the large N limit. The authors of ref. ||^ 
have calculated some correlation functions in two regimes: short range correlations (Aj — 
Xj) ~ 1/A^, and mesoscopic correlations (Aj — A-,) ^ 1/N. Note that the polynomials 
P/v and P/v_i oscillate at a frequency of order N (which corresponds to the discrete 
spectrum of a matrix of size N x N finite), and therefore, all the correlation functions 
will present such oscillations. 

In the short distance regime, these oscillations give the dominant behaviour, and 
eq. ([TTl|),([T|) lead to: 

sin [7rA^(A - At)p(A)] 



K(A,/i) 



7riV(A-/i) 



and all connected correlation functions follow from Eq. (| 1.3|) . All these features 
are studied in detail in ref. ^ , and we will now concentrate our attention on the 
mesoscopic case. 
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In the regime Aj — Aj ^ it is interesting to consider smoothed functions, 

defined by averaging the fast oscillations. For instance we find that: 

\n (\ n\] - ~^ ^ 1-C0S(^C0SV> 

L^2lA,/ij smooth — r.^T2 2 7\ ^2 ^ ^ 1 ' ^ -'-■-'-^j 

where 

a+6 h—a a+b b—a 
A = ^ ^cosy?, /^ = ^ ^cos?/-. (1.13) 

Brezin and Zee noted (ref. 0]) that the smoothed higher order n-points correlation 
functions vanished identically at the order l/N"^ for n > 2 (Note that k being of 
order 1/A^, the connected n-point function will be of order . Indeed, we will 

prove below by a direct method, i.e. without smoothing Eq. ( | 1.3|) , that they are of 
order 1/A^^'-"'"^''. The method is based on the study of the loop correlators, and these 
correlators have already been calculated recursively by the authors of [|l^ by the loop 
equation method. 

1.3 The n-Loop Correlation Functions 

Let us consider the following functions: 

u;„(^i,...,^„) =iV"-2/tr L- X ... X tr ( 1.14) 



Zi- M Zn- M 

They are related to the previous correlation functions by the relations: 

Un{z^, ...,Zn)= iV2"-2 / ]-[ _^ O^^X^^ ___^Xn) ( 1.15) 

ujn is called the n-loop correlation function, because it is the Laplace-transform of the 
partition function of a discrete random surface limited by n loops (the Zi are conjugated 
to the lengths of the loops, see p, ||, [1^] ). This remark allows to understand the 
topological origin of the factor iV^"^"^: indeed, the Laplace-transform of the complete 
n-point correlation function which is of order 1, would be the partition function of 
every surface (not necessary connected) with n boundaries. Each surface contributes 
with a topological weight where x is its Euler character. The leading term is 
the most disconnected one, with x = ^ (indeed, such a surface is made of n discs, 
each of them having x = 1), while the connected term has x = 2 — ri (it is a sphere 
(x = 2) from which n discs have been removed). Therefore, the relative contribution 
of the connected part to the complete n-loop function is ]\f^-^"-. As stated before, one 
can not only find the large N scaling of the connected correlation functions, but also 



compute them exactly by the recursive method of ||T0[. Let us give a brief summary of 
this method in the case which we are interested in. 

First remark that relation (| 1.15|) can be inverted: Un is analytical except when 
some of the Zi belong to the interval [a,b]. On can then be expressed in terms of the 
differences of Un between opposite sides of the cut. For instance: 

p(A) =^(a;(A + zO)-a;(A-zO)) 
02(A, fi) = -(2ih^{^2{'>^ iO, iO) -uj2{X + iO, /i - zO) - u;2(A - zO, yU + zO) 
+cu2(A-iO,/x-iO)). 
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and for general n: 



— i 1 

( 1-16) 

Note that these functions are directly the smoothed correlation functions, since we 
first compute Un in the large limit at complex arguments (which suppresses the 
oscillations), and then take the discontinuities along the cut. 

The main tool of |Tn] is the loop-insertion-operator, which acts on the free-energy 
of the Matrix-model, and gives the loop-correlation-functions: consider the partition 
function: 

-NiiV{M) 



Z = e'' ^ = J dM e 
with V{z) = J^kLigkz'' /k, and define the loop-insertion-operator: 



Note that with this definition 

bV{z') _ 1 _ 1 
SV{z) z — z' z 

The ujn are obtained from the free energy F by the repeated action of this operator: 

5 6 ^ 



1.181 



^n{Zi, . . . , Zn) 



6V{z^) 6V{z^) ' 

■^n~l\Zli . . . , Zn-lj- 



^V{Zn) 

(1.19) 

It is not necessary to calculate the free energy F, since we already know the one-loop 
function uj{z). We have already emphasized that oj{z) satisfies a linear equation Eq. ( 
1.1U| ), and thus, all its derivatives satisfy the same linear equation, with a different 
l.h.s. . Then, analyticity properties, and boundary conditions determine the form of 

The linear equation for uj is: 

cj(A + iO) +cj(A-zO) = l^'(A). (1.20) 
By a repeated action of the loop insertion operator, we obtain: 

UJ2{\ + IQ,Z)+UJ2{\-IQ,Z) = - . (1.21) 

[z-xy 

and for n > 2: 

uJn{X + iO, Z2,..., Zn) + a;„(A - iO, Z2, . . . , Zn) = . ( 1.22) 
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The function uj{z) has the form 

iviz) = \ {V'{z) - M{z)^{z-a){z-h)) 
where M{z) is a polynomial such that uj{z) ~ 1/z when z oo, and therefore: 

V'{z) 



M{z 

Then 



.\ji.z-a){z-h)^ 



p(A) = ^M(A)v/(A -a){h-\) = i^M(A)^ sin(^ , 

where we have used the parametrization ( | 1.13|) . 

The two-loop function is also completely determined by the linear equation ( 1.21|) 
and boundary conditions: 

1 1 / {x-yf -{x-a){x-h)-{y-a){y-h) 
^i^-yy\ ^{^^-a){x-h)^{y-a){y-h) 

and thus, in agreement with Eq. ( | 1.12|) : 

n (\ \= 1 1 l-cosy?cos^ 

^2KA,^l} 2Ar%2 _ sin(/?sinV^ ' 

The other loop functions all satisfy an homogeneous equation, and can be written: 

u;„(Ai, . . . , A„) = TT Wn{\i, A„), 

with now A„ = ^ — ^ cosy^n, and where the Wn are some symmetric polynomials 
of degree less than 2n — 5 in each Aj, which are no longer determined by the boundary 
conditions. It is necessary to directly use the recursion relation ( | 1.191 ). Since uj2 
depends on the potential V{M) only through a and b, we need the actions of loop- 
insertion operator on a and b, for instance 5a/6V{z), 6'^a/6V{z) 6V{z').... For this 
purpose, the authors of [1^ introduced the following moments of the potential: 

Mk = <pdz 



2m J (z — a)'' 



\J i.z - a){z - b) 



2m / {z- bY ^(^z^a){z-b) 

(the integration path turns clockwise around the cut [a, 6]). 
They are such that: 

M{z) = J2 Mk+i{z -af = Y. h+i{z - bf 

k k 
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The Mk and Jk are linearly related to the coefficients of the potential V, and if V is 
a polynomial of degree w, M is of degree v — 2, and there are only v — 1 independent 
coefficient among the and J^. Note also, that if V is even, we have = — (— 1)^ J^. 
The cut end-points a and b depend on V through the conditions that (see |10|, p^): 

^fd^ , ^'^-^ =0, 

y/iz-a){z-b) 

dz =2. 



2i-K 



y/iz-a){z-b) 

It follows (using Eq. ( 1.18|) and performing the contour integrals): 



Sa 
5V{z) 

Sb 
SV{z) 



1 1 



1 1 1 
•^1 ^-f> ^y{z~a){z-b) ' 



In order to determine the higher order derivatives of a and b, we need to differentiate 
the coefficients Mk and Jk- 



2k+l 5a 

5V(z) 2 SViz) V^^fe+l (z-a)'' ^ 

I 1 5b ( Ji Jj_ _ v^fc-1 Mi+i \ 

~^ 2 5V{z) \{b~a)l' (2-a)fe ^/=0 {b-a)k-l J 

and analogous formulae are obtained for the by the exchange a ^ b. 

With these tools, we can now determine the a;„ recursively. For instance is given 
in 13 : 



Ml 



u;3(x,?/,2;j 



a—b 



%{^^{x-a){x-b)^{y-a){y-b)^(z-a){z-b)^ 



Let us also write 074. For this purpose, we ffist define the polynomials 



Q{xi,a) = Yl 



[Xi 



i=l,4 



With this notation 



W2 



(-3(6- 



, M2 



m[Q{xi,a)Q{x„b)] 

+3-i.g(x„6) [(b-a) 



Q{xi, b) - 3(a - b)^Q{xi, a) 



+?>^Q{xi,a) 



[(xi — a){x2 — a){x'i — b){xA — 6) + 5 terms 



where the last additional terms symmetrize in the four variables. The connected func- 
tions On, are then obtained by Eq. ( | 1.16| ) and Eq. ( | 1.22| ), and they are simply given 
by: 



On{^li ■ ■ ■ , A„) 



-1 



— u;„(Ai + zO, . . . , A„ + zO). 
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1.4 Universality of the smoothed correlation functions 

- We observe that the only universal features of Oi{X) = p(A) are the square-root 
singularities at the edge of the distribution, otherwise the function p(A) is potential 
dependent. 

- The two point function is universal. 

- The n-point smoothed correlation functions with n > 3 can be calculated recur- 



sively by a systematic method [10|. The main property is that they consist in a sum 
of a finite number of universal functions and involve only the first n — 2 moments of 
the potential. That means for instance, that two potentials V and V* induce the same 
three-point function as soon as they yield the same Mi and Ji coefficients, but they 
don't need to be identical. However, the n-point correlation function is no longer given 
by a gaussian model, since it is of order l/A^^""^. Perturbative corrections to the gaus- 
sian model have to be considered. Note that the determination of correlation functions 

allows an evaluation of the moments Mk, in a case where the potential V is unknown. 
Let us recall |1T0| , [Tl[| that these moments play an important role in the study of critical 



points. Since 

p(A) = LMiX)^{X-a)ib-X) 

we see that if some moments vanish, the behaviour of p near the end-points is no 

longer a square- root. When the m first vanish, one finds p ~ 

corresponds to the m*^ multicritical point of the one-matrix model of 2D gravity H, |T| . 



1.5 Conclusions 

In this section we have recovered, by a completely different method, the results of 
ref. 0] concerning the two-point eigenvalue correlation function of a random hermitian 
matrix in the limit in which the size N of the matrix becomes large. This method of 
orthogonal polynomials will be extended to another model in the next section. However, 
this method is not convenient to evaluate smoothed correlation functions. The way to 



do that, is the recursive method of [10|. Note that this method can easily be extended to 



the 0{n) matrix model using the two-point function and the genus expansion procedure 



presented in [|^, [T^ . 



Brezin and Zee in ref. 0] have shown that the two-point function is universal, and 
therefore identical to the function of the gaussian matrix model. The gaussian model 
is the fixed point of a renormalization group ^ and a direct RG analysis should be 
performed to put this result in perspective. 

In the same way, Brezin and Zee have shown that higher smoothed correlation 
functions vanish at leading order. The contributions calculated by the loop insertion 
method of [|10[ should be considered as corrections to the leading scaling behaviour. 



The explicit expressions show that they depend now on successive moments of the 
potential, indicating an implicit classification of the deviations from the gaussian model 
in terms of their irrelevance for large. Here also it would be interesting to confirm 
the qualitative aspects of these results by a direct RG analysis. 
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2 The 2-matrix model 



The previous method to evaluate the correlation functions of the eigenvalues of a single 
random matrix can be extended to the case of a multi-matrix model. It is known from 
2D quantum gravity that the multi-matrix model provides a larger class of critical 



points, and therefore, the universality properties may be different from the one-matrix 
case. It has also been established |[T^ that a two-matrix model is sufficient to reproduce 
all the critical behaviours of multi-matrix model, we shall therefore concentrate our 
attention on this case, but all the method may be generalized to any chain of matrices, 
provided you can use the Itzykson-Zuber integration to define orthogonal polynomials 

A multi-matrix model can be regarded as a particular one-matrix model with a non- 
polynomial potential, if you integrate over all the other matrices. In this case, you are 
interested in the correlation functions of the eigenvalues of only one of the matrices, and 
we will see below that in the large limit, we recover the same universal behaviour as 
for one matrix with polynomial potential. The multi-matrix model presents also some 



interest in the study of coupled systems |^ |T8|. In that case, you are interested in the 
correlations of eigenvalues of different matrices. We will see, that in the large N limit, 
the eigenvalues of the two matrices are essentially uncorrelated. 

First, let us define the model and the problem. We will try to use the same presen- 
tation and notations as for one matrix: 

we consider two stochastic matrices M and M hermitian, of size N x N, with a 
probability: 

1 -Nil (V(M) + V(M) - MM) , ^ 

V{M,M) = -e ^ ^ ' ^ ' ' (2.1) 

Where V and V are two given polynomial potentials, not necessarily of the same degree. 
This probability distribution induces by integration over the angular degrees of freedom 
of the matrices (the unitary group) a joint probability density for the eigenvalues of 
both matrices. Using the famous formula of [19| , the result reads: 



P,A..A,) ^ iA,A)A(A)det,e^^-^^)ne-^<^'<''' + ( 2.2) 

^ i 

where the A are the Vandermonde determinants: 

A(A) = n(A. - A,) , A(A) = n(A. - A,) 

i<j i<j 

The quantities we are interrested in are the correlation functions of eigenvalues. We 
consider the operators: 

0{X) =itr5(A-M) =±j:^6{\-X,) 

OCX) =^tr5(A-M) =^^^5(A-A,) ^ 

One would like to evaluate mean values as: 

Oi(A) = (0(A)) 
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' conn 

m 



OliX) = {OiX)] 
02,o(A,/x) = (0(A)0(^)), 

n0(Ai)n^M) • (2.4) 

i=l 1=1 / conn 

These correlations are related to the partially integrated eigenvalue distributions: 

Pn,m(Al, . . . , A„, Ai, . . . , Xm) = J ■ ■ ■ J Yl '^i Yl ^(Ai, . . . , Aat, Ai , . . . , Xn) ( 2.5) 

i>n j>'m 

Indeed, when all the Ai, . . . , A„ and the Ai, . . . , A^ are distincts, we have: 

N\ N\ 



n,m 



Unlike the one-matrix case, we have here three kinds of correlation functions: those 
with only eigenvalues of M, those with only M, and the mixed ones. For instance, the 
densities of eigenvalues: 

p(A) and p(A) , 

and the two point-functions: 

P2,o(A,/i) po,2(A,/i) pi,i(A, A) 

From the one-matrix model case, we guess that the function p2,o(A,p) shows an inter- 
esting behaviour at small distance A — p ~ but what about the mixed correlation 
? A and A are not quantities of the same nature (A scales as 1/A), and nothing will 
happen when A is close to A (except perhaps in the symetric case V = V). Also, 
as in the one-matrix case, another interresting problem is to calculate the smoothed 
correlations at large distances ( ^ ^/N). 

As in the one-matrix case, the method of orthogonal polynomials provides the short 
distance correlations in terms of some kernels (here there are four kernels), and these 
expressions can be smoothed to obtain the smoothed two-point functions. The higher 
order smoothed correlation functions can't be obtained by this method, but by the 
loop equations, which are too complicated in the two-matrix-model, i.e. there is not 



yet any generalization of |10 



First, we recall the method of orthogonal polynomials and introduce the kernels. 



2.1 orthogonal polynomials 

We will consider two families of polynomials P„ and -P„ (n is the degree) with the same 
leading term, and which satisfy: 

HXHXPa^P (V, -N{V{X)+VCX)-XX) 

dXdX Pn{X)Pm{X)e = On,m (2.6) 
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We also define the wave functions: 

V^„(A) = P„(A)e~^^^^^ 

Wliicli we note: 

< n| = ^n(A)^< A| 
|n > =ip(\)\\> 

tliey form two dual spaces with the bilinear form: 

< n\fn >= [ dXdA^„(A)?/im(A)e^'^'^ = 



2.7) 



2.81 



with this duality, it is possible to map each of these spaces into the other one, we thus 
define two other families of functions (which are not polynomials): 

Xn(A) = /dAVi„(A)e^^^ , . 

A^AA ^ ' 

Xn(A) =/dX?/^„(A)e 

As in the one matrix-case, partial integrations over eigenvalues in Eq. ( | 2.2| ) are per- 
formed by expressing the Vandermonde determinants in terms of orthogonal polyno- 
mials, and result into combinations of four kernels: 



( 2.10) 



We then have (see Appendix 4 or E^]) 



p(A)=if(A,A) p(A) = i7(A,A) (2.11) 



pS:!(A,A) =-i^(A,A)K(A,A) + ^ir(A,A) 



2.12) 



As for one matrix, we have to evaluate these kernels, this can be done by a generalisation 
of the Darboux-Christoffel theorem. We have for instance: 



("41) 



—NXX if ~ degy'fc-l 

K{X, A)e = JTt \ '^i'N'ipN - i'N-k+i'ipN+i I ( 2.13) 

^ k=l i=0 



wich involves only ijjn and ipn with n—N ~ 1. Note that this is a differential equation for 
K, but in the large A^ limit, the partial derivative will be replaced by a multiplication 
by a function A (A), as we will see below. Such Darboux-Christoffel like formulaes 
exist also for the other kernels. Therefore, our next task will be to find an asymptotic 
expression for the ipn in the limit A^ large and n — A^ ~ 1. 
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Let us now explain the method (similar to the one used to calculate the two-loop 
functions). We consider the following operators of our Hilbert spaces: 

f-lc? ~ ~ 1 d , ^ 

d . -6 



such that: 



n ,(p = -f- , X 
dn 



< n\\ = XipnW < A| 

X\h > = Xipn{X)\X > 
n\m >= m\m > , x\n >= |n — 1 > 
< mlfi = m < m\ , < n\x =< n + 1\ 



It is easy to see from the degrees of the poljTiomials that A and A have the form (for 
instance, a multiplication by A can raise the degree of at most 1): 



2.15) 



2.17) 



A = a{n)x + J2k>o(^k{n)x~'' 

the equations of motion ( just integrate Eq. ( | 2.6| ) by parts) give simply: 

P = -X , P = -X (2.16) 

and: 

p =-nA) + ^xf +o(x2) 

The expansion of Eq. ( | 2.17| ) in power series of x, using commutation relations like 

[x, a(n)] = {a{h + 1) — a{n))x 

gives a set of recursion relations between the ak-,ak- these last equations allow (in 
principle) to determine all the ak and a^, by induction. Note that only a finite number 
of ak are non zero, and more precisely: 

degV" 

X = a{n)x+ ^ ak{n)x~^ 

k=Q 

^ degV 

X = x^^a{n) + ^ x^ak{n) 

k=0 

Another remark is that all these recursion equations are equivalent to the cannonical 
commutation relations: 

M = ^ (2-18) 

Up to now we have written only exact equations, let us now discuss the large limit 
approximation . 
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2.2 Large N limit 

It is the classical limit, i.e. we can drop the commutators, and replace operators by 
numbers. We will also be interested in the limit n — ~ 1, and we will most of time 
replace the operator fi by the value A^. The ak{n) are then numbers ak which no more 
depend on n. We thus write: 

A = ax + Efcio "fcX " ( o ^Q\ 

A = ax + 2Zk=o 

At leading order the recursion equations for the become algebraic equations, they 
come from the expansion in power of x of the classical equations of motion: 



V'{\)-\ =xl + 0{x 



( 2-20) 



still, we remark that these algebraic equations can also be obtained from a Poisson 
bracket: 

(dXdX _ dXOy 
\ dn dx dx dn 



;)={A.A} = 1 (2.21) 



The left member is expanded in powers of x. In particular the coefficient of x° is: 



ia^ — kakCHk) = — tt 



which gives after integration: 

a^ -Y.kakak = -- = -l (2.22) 

From these algebraic equations, one determines all the coefficients a^, and obtains two 
functions \{x) and A(x) of an auxiliary variable x\ 

A(x) = ax + Y.k>o(^kX~'' ('2 23) 

A(x) = ax^^ + Yjk>QOikX^ 

In particular, we can define a function A(A): we fisrt get x from the first relation and 
insert it into the second one. Such a function is multivaluate, it has a cut along a 
segment [a, h] of the complex plane. Actually, it can be seen from the loop equations 
(see appendix 6 or [^] ) that V^'(A) — A(A) = a;(A) is precisely the resolvent: 



We now arrive to the aim of those calculations, i.e. finding an asymptotic approxi- 
mation of the orthogonal polynomials, and kernels. 

The leading behaviour of ^/'n(A) when A^ ^ 1 (classical limit) and n ~ A^, is simply 
given by: 
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indeed P = —A according to Eq. ( | 2.161 ) , and therefore: 



dX 

As in the one matrix case, we define: 



-A^A(A) 



((A) =i J A(A)dX ( 2.25) 

and obtain: 

ipn{X) ~ e 

This rough approximation is not sufficient to determine the kernels, we have to go to 
order 0{N^). Remark that the dependance in — n is easily found, since: 

< n\ =< Ar|x"-^ 

or in other words: 

A proof is presented in Appendix 5, we give here only the final result, which is very 
close to the one obtained for one matrix (remember Eq. ( | l.llj ) ): 

^n(A) JiEx -jw^ ^ 



dr 



PT 1 -^/r A(z/)A'(y)ci/ ^ 

— 

2^ 1 N j:^X{y)X'{y)<ij 



2.26) 



XniX) ~y'-Ex^=^e ° X 

V dc 



. 1 Nj:^X{y)X'{y)dj 



dc 



Remarks: - The means that you have to sum over all the values of x which give the 
same value of A = A(x) (in the one- matrix case, we had only two values of x complex 

conjugate of each other (x = e"*"*"^), that's why the sum over x was replaced by a real 
part). 

- ip can be obtained from ip by the exchange x ^ 1/x, and ak ^ otk- 

- The expressions of x and x are obtained by a steepest descent approximation in 
Eq. ([g . 

- The lower bound of integration xq can be chosen arbitrary but must be the same 
for all. It can be chosen for instance such that A'(xo) = 0, and such that X{xq) = a is 
one of the cut end-points of the cut [a, h] of the resolvent uj{z). 
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2.3 Darboux-Christoffel theorem for the kernels 

We have already asserted, that there exists a kind of Darboux-Christoffel theorem, 
which gives the kernel (| 2.10|) in terms of a litle number of ipn- For K and K, we have 



only differential equations, but in the classical limit, the derivatives disappear. For 
instance, we have: 

= jfj:n=i r\N>®<N\x-^ 

with the operators x and y acting on the < N\ and \N > respectively, they commute, 
and we can sum up the geometrical serie: 

K{X, X)e'^^^ = ^j^\N X N\ ( 2.27) 

y 

in the classical limit, x and y become x = x(A) and y = x{X). Let us multiply this 
expression by A — A (A) which is: 

/ k-l \ 



X - A(A) = A(a;) - A(y) = (x/y - l) Uy -J^a^Y^ x^-%-' 

\ k j=0 



Therefore: 



ir(A,A)e — «y-E«.E^^'VM|iV><iV| (2.28) 

which is formally identical to Eq. ( | 2.13| ). By Laplace transforms, we can express the 
other kernels: 

11/ \ 

H{\ /^) = J^Y^ \ tpN-k+iWxN+i{f^) 1 ( 2.29) 

^ \ k 1=0 / 

and analogous expressions for H and K. 

We can now insert ( | 2.26] ) in ( |2.28| )... ( | 2.29] ), and compute the correlation functions 



m, ([HI)- 



2.4 Correlation functions 

Let us first compute the kernel H{X,fi), and in particular the limit X = fi which will 
allow to fix the normalizations. From ( | 2.27| ) we have: 



H{X{x),fi{y)) = _i ^n{X)xn{^) 



and with ([2:261) 
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All the X and y such that A = A(x) and yU = \{y) contribute to the result in the general 
case. But in the case X — fi small, the terms such that x — y is small are dominant , 
and there are two such terms, complex conjugate of each other. We thus have: 

^ R ^ 1 -iVA(A)(A-/i) 

Since we have uj{\) = V'{\) — A(A), 

-lmA(A) = Imcj(A) = -7rp(A) 

we recover if (A, A) = p(A). (Actually, this relation was used to fix the prefactor 
cte= 1/^2^ in (|T2|) ). 

The short distance (A — yU small) expression of H is identical to the one obtained 
in the 1-matrix case: 

, sin7rAr(A-/x)p(A) -N{X - fi)Re~X{X) 
vrAr(A-/.) ' 

The two point function can then be derived: 

pg(A,/i) = -/7(A,/i)if(/x,A) = V(A) (^)' ( 2.30) 

with e = A^(A— /x)p(A). This result confirms the universality of the two point correlation 
function in a case more general than 0] . From this expression of H and from Appendix 
4, we can also find all the connected correlation functions with more than two points, 
in the short distance regime. The result is still identical to the one-matrix-model's one. 

Now let us come to the new aspect of this model: the mixed two point function. 
We will see that it is always zero at the order we consider. We have from ( | 2.27| ) and 



2:261) : 



K{X,X)= V e 

x,x',y,y* ^TtN X-y J_^S^ 

dr dy 

K is proportionnal to and vanishes at large N, except when x —>■ y. The corre- 
lation function thus vanishes as l/N"^. 

The only domain in which this correlation could be larger stands along the curve 
A = A(A) in the (A, A) plane (i.e. when x = y). Unfortunately, this curve lies outside 
the interesting region A G [a, 6] and A G [a, 6]. So, we find that in the large N limit, 
the two matrices are uncorrelated. 

However, let us suppose that: x — y <^ 1. Then: 

1 -N^^X'{x)X'{x) 

K{X{x),X{y)) , . e 

2nN{x -y)^~X'{x)X'{x) 
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if we set e = \/N (x — y) \J —X'{x)\'{x) we obtain: 



1 1 eV2 
~ Re — ^-e 




A-A(A)= - 

pS(a,a) ^ 



We still observe some kind of universality, but the meaning is not clear, since it does 
not concern the true physical eigenvalues. 

2.5 Smoothed two point functions 

When A — is not small, as for one matrix, the correlation functions are not universal. 
However, we see that equations (| 2.26| ) contain very fast oscillations of frequency N . 



At large distances, these oscillations can't be observed, the physical quantities are the 
smoothed functions. So, inserting expressions ( | 2.26|) in (| 2.101 ) and then in (| 2.12|) 



, we find the complete two point correlation functions in the large N limit. Then 
we smoothing the large frequency oscillations, i.e. suppressing the terms of the form 
exp / AA', we get: 

P(2!o)(A)p) = -:;^2ff2^Y.\[x)=\,\{y)=f^ \'(x)\'{y) (x-yY 

P{0,2)(A,/i) = ;^27v2 4Ea(x)=A,A{j/)=A A'(x)A'(y) (2-31) 
P(1,1)(A'A) = :;^^Y.X(x)=\:\(y)=\ X'{x)\'{y)'(FW 

As in section (1.3), such smoothed correlations could also be obtained by the loop 



correlators. The authors of JT^ have calculated explicitely the two-loop correlators 
and found: 

c^o,2(5i,52) =-g,,g,,ln -(^:y) (2.32) 



ZI-Z2 

Where the functions x{z) and x{z) are defined by inverting (| 2.23|) . Both results 
coincide. 

As for one matrix, the other smoothed correlation functions can't be obtained by 
this method. Indeed, since the kernels K, H, K, H are of order 1/A^, we can compute 
the n-point function only at order while we have seen by a topological argument 

that the n-loop function is of order l/A^^""^. To continue, it would be necessary to 



generalize the methods of to the multi-matrix model. However, this has not been 
done yet, but we can guess that the higher order smoothed correlation functions can 
be expressed in terms of the functions x(A), and an increasing number of coefficients 
depending on the potentials V and V. They should present more universality than the 
1-point function. 
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2.6 Conclusions 



By the analysis and the exact calculation of the correlation functions of the 2-matrix 
model, we have confirmed with more general hypothesis the universality properties 
that Brezin and Zee have established for a 1-matrix model. When one considers only 
one of the matrices, every thing seems to happen as in the 1-matrix case: the short 
distance 2-point correlation function (| 2.30|) is exactly the same. But the long distance 
two point function is no more universal (cf [^]): it can be expressed in terms of an 



auxilary function x(A), but in the general case, it is not as simple as (| 1.12|) . The 
functions x{X) and x{X) themselves, are not universal, they depend on the coefficient 
a, ak, oik (one of them can be eliminated by Eq. (| 2.22|) ), i.e. degl^ + degV^ coefficients. 
In other words we have just performed a convenient change of variables among the 
parameters gk,gk- 

The new aspect of the two matrix model, is that we can check the correlation 
between two coupled random matrices. We then find that in the large limit, they are 
uncorrelated. Indeed, the correlation function is most of time of order l/N"^. The only 
domain in which this correlation could be larger is not physical, i.e. when A and A verify 
A = A(A) (we have seen that if A G [a, 6], A(A) is not real, indeed ImA(A) = 7rp(A)). 

We have not treated the case of a chain of random matrices coupled by a weight of 

the form e ^ , but it seems reasonable to think that all the method developed 

in this article should apply with no difficulties to it , and it will be done in a following 
paper [^. We would then find the same general feature: the short distance correlation 
function of one of the matrices is still (| 2.301 ) . The smoothed functions still look like 
( I 2.31| ) , but with more complicated functions x(A), involving a set of coefficient 
equivalent to the set of the coefficient of the potentials. 
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Appendix 1. 

Orthogonal polynomials P^: an explicit expression. 

Let us show that the orthogonal polynomials Pn defined by the orthogonality con- 
dition: 

r -NViX) 

{Pn ■Pm)=Snm = J P„(A)P„(A). 

are given, up to a normalization, by equation ( | l.Tp : 

Pn{\)(x (TMe ^ ^det(A-M). 

First, this integral clearly yields a polynomial of degree n in A. Let us then verify 
the orthogonality property: after integration over the unitary group, the integration 
measure / dM reduces to an integration over the eigenvalues of M, and the Jacobian 
of this transformation is a square Vandermonde determinant: 



^n(A)oc/ n dXie ^ '^A-A,) A2(Ai,...,A„) 

•' i=l,...,n 

where 



A(Ai,...A„) =n(A.-A,). 



In this form we recognize a more classical expression p5[. Then, setting A = Aq: 

-NV{\ 



< Pn{Xo) ■ A^ >oc 1 II dX,e ^ A(Ao, Ai, . . . , A„) A(Ai, . . . , A„)A^ 

i=0,...,n 



The first Vandermonde is completely antisymmetric in the n + 1 variables, we can 
therefore antisymmetrize the factor A(Ai, . . . , A„)AJ}^, the result is zero if m < n because 
the only polynomial completely antisymmetric and of degree less than n — 1 in Aq is 
zero. Thus 

< P„ • A" >= . 

it proves that P„ is orthogonal to any P^ with m < n, and that < Pn-Pm >= as 
soon as m ^ n. 

Remark that a similar integral representation can also be found for multi-matrix 
models. Consider two families of orthogonal polynomials P„ and P„, such that: 

< P„(A) ■ P„(A) >= I jdX dA e~^^^^^^ + ^^^^ ~ ''^^^P„(A)P^(A) oc 
Then, we have: 

f f ~ -N tr \V(M) +V(M) - cMM] 
Pn{X)= (MdMe L I I J Met(A-M) 



where M and M are hermitian n x n matrices. Similarly: 

~ ~ f f ~ -N ti [V(M) + V(M) -cMM] 
P^(A) = / (M (M e ^ ^ ^ ^ ^ Met(A-M). 
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Appendix 2. 

Connected correlation functions and the kernel «:(A,/i). 

Some exact expressions exist for the correlation functions of matrix models [T^j, in 
terms of a kernel K(A,/i). We again consider the matrix distribution ( | 1.1|) : 

, , 1 -A^trV^(M) 
P(M) = -e ^ ^ 

The corresponding measure can be rewritten in terms of the eigenvalues Aj of M and 
a unitary transformation U which diagonalizes M: 

i=\...N 

(A = Y[i<^j{\i — Xj) being the Vandermonde determinant). Therefore, the probability 
that the eigenvalues of M are Ai, . . . , Atv is: 

PAr(Ai, . . . , Aat) oc Yl 6 A^(Ai, . . . , Aat), 

i=l...N 

and the correlation functions are obtained by partially integrating over some eigenval- 
ues: 

P(Al) = / Ili=2,...,N dX i pAf (Ai, . . . , Aat), 
P2(Ai, A2) = / ni=3,...,Af dX i pAf (Ai, . . . , Aat) 

... and so on. 

The Vandermonde determinant A can be written as: 

A(A,) = det Ar\ 
and thus, after some linear combinations of columns of the matrix: 
A = det nj_i(Ai), n„(A) = A" + O (A""^) , 

identity true for any set of polynomials {H^} normalized as above. In order to perform 
the A integrations, we choose n„ oc P„, P„ being the orthogonal polynomials (| 1.6|) . 
A^ is the product of two such determinants, therefore it is the determinant of a matrix 
product: 

A^ (X det K{Xi,Xj) 

Kix„x,) =j:^:=o'Pkix.)Pkix,). 

(2.33) 

The proportionality constant in Eq. ( 2.33| ) is here irrelevant because the eigenvalue 
distribution is normalized. The Darboux-Christoffel formula (appendix 3) tells that: 

. ^ ^ PN{X)P^^^{^l) - PMPn-i{X) 
A — p 
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where a is a normalization constant depending on A^, and a — {a — h)/A when N — > oo. 
The important properties of K{\ fi) are: 

e'^^^""^ K{\,u) =1, /d/ e~^^^''^i^(z/,z/) = iV, 
-NViv) 

/di/e ^ 'k{\u)K{p,h) =K{\h). 



An exphcit expression for p„ 

-NV{Xi 

i=n+l,...,N i=l,...,N 



Pn(Ai, . . . , A„) = -^y" n n ^ ^^^^^^ detK{Xi,Xj), 



can be obtained by successively integrating over eigenvalues. Using the rules of K 
integration it is easy to prove by induction 

7V"(A^ -n)! 
Pn (Ai, A2, . . . , A„) = — det K {Xi, Xj) , 

where we have introduced the reduced function 

1 -iN/2)[V{X) + V{,,)] 
/«(A,/^) = -^e K{X,ii). 

Therefore, when the Aj are all distinct we have: 

(C»(Ai) . . . 0{Xn)) = det k{K, a,) = Y.i-^y n 

<T i 

The connected function will involve only the sum over the permutations a such that 
Yii ii{Xi, Ao-J cannot be split in the product of two cyclic products, i.e. only cyclic per- 
mutations will contribute to the connected function. This intuitive result is a classical 
combinatorial identity, can be proven fromThe connected function can thus be written: 

(0(Ai) . . . 0{Xn))conn = (~l)"^^n ^permutations a Hi l^iKi, Ki+i) , 

= (-1)"+! [k (Ai, A2) K (A2, A3) ... « (A„, Ai) + • • •] , 

where the additional terms in the r.h.s. symmetrize the expression. Or in a more 
compact way: 

(o(Ao . . . o(A„))_ = (-1)'^-^^ y: n i>^i^ K) 

cyclic permutations a i=l 

Appendix 3. 

Derivation of the Darboux— ChristofFel formula. 

The polynomial AP„(A) can be expanded on the basis of the Pm with m <n + 1: 

n+l 

A-Pn(A) = ^ QnmPm (A). 

m=n— 1 
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The orthogonality condition ([ 



1.6| ) imphes that the matrix Q is symmetric: 



Qnm ^ ' Pm ^ ^ Pn ' ^Pm ^ Qmn • 

The polynomial (A — fi)K{X,fi) can thus be written: 

N-l / N N \ 

k=0 \i=0 j=0 J 

All the terms cancel, except the upper-bounds: 

(A - /i)ir(A,/i) = Q^,^_i(P,v(A)P^-i(/i) - P7v(/^)P^-i(A)), 

therefore a = Qn,n-i- In the large limit, it is possible to calculate a. The simplest 
way of doing this is to calculate APjv„i from expression ( | 1-1 Since 

f a + b b — a \ a + b b — a 

I — — cos ip I cos [ip — ip) = — - — cos [ip — ip) — (cos ip + cos [ip — Zip)) , 

we have: ^ ^ 

APat-i = — - — Pn-1 1 — {Pn + -P/v-2) 

and therefore a = {a — b)/A. 
Appendix 4. 

Partially integrated distributions of the 2-matrix-model. 



Exact expressions for the correlation functions have been recently derived in 
. Hereafter, we just give the results and present a diagrammatical method to express 
the partially integrated distributions of type ( | 2.5| ) in terms of the kernels K, H, K. 
For instance, generalizing the one-matrix result, we have: 

1 " 

-p„,o('=)(Ai, . . . , A„) = (-1)"- E n H{K. 

^ cr i=l 

We observe that p^^^^ is a symetric sum of cyclic products of kernels H. The general 
pI^^j^ have the same form, except that the other kernels appear. The method consists 
in writting all the possible permutations of the eigenvalues Ai, . . . , A„, Ai, . . . , A^, and 
multiply the kernels of two consecutive eigenvalues in a cyclic product. All this can be 
represented in a diagramatic method: 

-First write a permutation of all the n + m eigenvalues as a closed oriented chain. 

-assign to each link between two consecutive eigenvalues a factor according to the 
following rules: 

Aj — > \j = —H(\i,\j) ^ Aj — > \j = —K(\i,\j) 

Aj — > Xj = —H{\i, Aj) Aj — > Xj = jv — K{Xi, Xj) 
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-then — p'-'^-' is the sum for all these diagramms of the products of the links. You can 
verify the number of terms in the sum is 

min(n,m) i i 

^ nl ml 

For instance we have already given the one and two point functions in ( | 2.12| ), we 
write also here the 3-point functions: 

pgi)(Ai, A2, A) = kCX, Xi)H{X,, X2)K{X,, A) + KiX, X2)H{X2, Xi)K{X2, A) 
-j^H{Xu X2)K{X2, A) - j^H{X2, Xi)K{Xu A) 



Appendix 5. 

An approximation for the two matrices orthogonal polynomi- 
als. 

The proof follows the same method as in the one-matrix case: we have the integral 
representation (Appendix 1) 



. -A^tr (V{M) + V{M) - MM) 

Pn{X) = cte / dM"^"dM"^"e ^ ^ det(A - M) 

wich we write: 

P„(A) = Z(^A,^ = — ,/i = — j=e 

tr fV(M) + V{M) - mm) 
Z(A,^,/i) = / dM"^"dM"^"e ^ ^ ^ 

V{z) = V{z) -h\nX-z 

We observe that 

dF{X,g,h) \ \ h\ 

n — = u]{z = A, A, g, h) 

dX 

where oj{z) is the resolvent corresponding to Z{X, g, h). We need the leading behaviour 
as A^ — ^ 00 and A^ — n ~ 1 of n'^F, and we need only its dependance in A. 

n- = Ar,c. 

We have to expand u in the limit g 1 and h —>■ 0: 

gu = Uo + {g- l)ng + hQh + O(^) 

with 

g = l,h=0 



^ _ dgLV 

uoiz) =uj{z,g = l,h = 0), o _ 

^''h — dh 



g=l,h=0 
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Then: 

dF 

n^— = NtooiX) -{N- n)ng + + 0{l/N^) ( 2.34) 

CIA 

We have aheady noted that: 

Uo{z) = V'{z) - ~X{z) 

and we can determine flg^h by the remark that a variation of the potential by an amount 
of 5V{z), induces a variation of u given by the two-point functions 002 = duj/dV: 

6lj{zi) = / ^^2,0(2:1, Z2)SV{Z2) + UJi^l{Zi, Z2)6V{z2) ( 2.35) 

Qg corresponds to the variation of the potential 6V = 1/gV — V ^ —{g — 1)V, and 
Qh to 6V = —hlnX — z. Inserting the explicit explicit expressions ( | 2.32|) of [|l^ into 
( I 2.35[ ) , and integrating by parts and taking the residues, we obtain: 

at 2; = A we have: 

Id dxjX) 

Remark that for Qg there is a simpler method: we have ipn = x"'~^'4'ni the (A^ — n) 
dependance is thus x""^. Therefore, we know that 

l],(A) = ^lna;(A) 

Inserting these expressions in ( | 2.34|) and integrating with respect to A, we obtain 
the asymptotic expression of ipni.^)'- 

fd^ -N r X(X)dX 
V'n(A) ~ cteW^^e ^ ' x"-^ ( 2.36) 

V dX 

Unfortunately, this result is obviously wrong. Let us remember that in the one-matrix 
case, at this point, we had considered the real part of this expression, and replaced 
the exponential by a cos . Another way to understand what happened, is that in the 
classical limit, we have replaced operators by numbers. A more carefull analysis would 
have given instead of ( | 2.24|) a differential linear equation of order degV^ for ip^) and 
thus we should add independent solutions which all have the form (| 2.36|) , but with 
the different values of x satisfying Eq. ( 2.23|) (Thanks to F. David for this remark). 



^n(A) ~ 2^ ctcfc W — e ( 2.37) 

The function x(A) is defined as the physical one, i.e. such that V'{X) — A(x) = uj{X) 
is the resolvent and behaves as 1/z for z — > oo. It is the one such that A(x) ~x'^oo oix. 
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At large A, only this solution contributes to ( | 2 .361 ) , while along the cut A G [a,b], all 
the X contribute. We then have: 

V'n(A) ~ cte2 2^ W— e x . (2.38) 

The same method gives also ipn, and to obtain the Xn you just integrate ( | 2.9| ) by 
steepest descent. Note that the cte factor can be fixed by the normalizations. At the 
end, you find i\22§^) . 

Appendix 6. 

The loop equations for the two matrix model. 

The loop equation refiect the invariance by a change of variable into the partition 
function: 

First consider the change of variable: 

M^M + e- ^ ^ 



'z-Mz-M 

At leading order in (i.e. approximating < TrTr > by < Tr >< Tr >) we get : 

(y'{z) - Lj{z) - Z) LJ{Z, Z) = ^PijZ'uJj{z) - LU{Z) 

where we have defined 

uj{z,z) = ( Tr \^ ^ . ) , uJz) = (ti^—mA 
^ ' ^ \ z-Mz-M/ n ^ \ I 

'^"^^'-^ V'(x\ - V'{v\ 

and P{x,y)= E Vr^^W = -^^^-^ 

Making the symetrical change of variable for M we find that: 

^ J2ijPijZ'uJj{z) ~Uj{z) ^ J2i,jPijZ'uJj{z) -Ljjz) 

V'iz)-iuiz)-z ~ V'{z)-Cj{z)-z 
then reducing to the same denominator, we get that 

{y\z) — ijj{z) — z^ ^^pijZ^ujj{z) — zuj{z) + V'{z)uj{z) = idem with z ^ z 

The left hand side is a polynomial in z of degree degl^' — 1, while the right hand side 
is a polynomial in z of degree degy — 1. Therefore Both sides are polynomials, let 
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us call this polynomial Q{z, z). This allows to write some algebraic equation for uj{z) 
Moreover, we have an expression of u{z,z): 

. . , Qiz,z)-{V'{z)-z){V'{z)-z) 

Uj[Z, Z j — 



{V'{z) - uj{z) - z){V'{z) - - z) 
The poles must cancel, so: 

Q{z,\{z)) = {y\z)-\{z)){y'{\{z))-z) when z = \{z) = V {z) - u:{z) 

Q{\{~z),~z) = {y\X{~z))-~z){y\~z)-\{~z)) when z = \{~z) = V' {~z) - Cj{~z) 

These algebraic equations can be shown to be equivalent to the equations of motion 
2 .201 ) found by the method of orthogonal polynomials. 
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